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Abstract. We report on investigations of the applicability of non-relativistic constituent quark models to the low-energy 
nucleon-nucleon (NN) interaction. The major innovations of a resulting NN potential are the use of the 3 Po decay model 
and quark model wave functions to derive nucleon-nucleon-meson form-factors, and the use of a colored spin-spin contact 
hyperfine interaction to model the repulsive core rather than the phenomenological treatment common in other NN potentials. 
We present the results of the model for experimental free NN scattering phase shifts, S-wave scattering lengths and effective 
ranges and deuteron properties. Plans for future study are discussed. 
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INTRODUCTION 

The NN interaction has been the object of study for over 50 years. The development of the highly accurate potentials in 
1990s showed that the NN interaction can be quantitatively expressed very precisely although the physical mechanism 
of the interaction has yet to be understood. Most of these models use a partial-wave decomposition of the NN 
potential and vary potential forms and parameterizations in dependence on partial wave channel. The NN potentials 
are either thought to rise from one-boson exchange (OBE) (e.g. Nijmegen and Bonn families) or are constructed of 
different functional forms multiplied by general operators (central, tensor, spin-orbit etc) which are then parameterized 
(Argonne and Urbana families), with one-pion exchange included to model the long-range part of the potential. Chiral 
perturbation theory is a modern alternative approach which attempts, to derive the NN-force from contact interactions 
and pion exchange as the low-energy QCD degrees of freedom. The Nijm I[l], Argonne vl8[2], N 3 LO[3], and CD- 
Bonn[4] potentials are used in this work for comparison. 

In this work, we present a potential motivated from the Non-Relativistic Quark Model, the Oxford potential. It is 
shown to be able to well reproduce NN scattering phase shifts data and the properties of the deuteron. 

This work is organized as follows. First, we briefly review two quark models of strong interaction dynamics. 
Secondly, we create the Oxford potential using these models. Then the results of our parameter tuning are presented 
for NN data and we show the high momentum behavior of the potential. Finally, we give our conclusions and plans 
for future work. 



QUARK MODELS OF STRONG INTERACTIONS 
3 P Model 

It was first suggested by Micu[5] that hadron decays might be dominated by the production of a qq pair from 
the vacuum with J PC = ++ quantum numbers, i.e. they are in a 2s+1 Lj = 3 Po state. Her idea has been applied to 
successfully estimate 3-point vertices with one additional meson in the final state[6]. 
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FIGURE 1. A diagram illustrating the technique for calculating quark model form-factors for use in meson exchange potentials. 



FIGURE 2. A representative diagram for the quark model OGE potential. The x denote the colored spin-spin contact interactions. 



Fig. 1 illustrates how to parameterize the effective field theory from the microscopic decay diagram using 
the 3 Po model. External hadron simple harmonic oscillator SHO wave functions are attached to the microscopic decay 
diagrams to account for non-perturbative hadronization effects and yields the following expressions for the form-factor 

J?' 

a, j3 are the SHO wavefunction widths for the nucleons and meson respectively and pi, pf are the initial and final state 
momenta. 



One Gluon Exchange 

Attempts to model the NN potential via a OGE interaction between bound state hadrons dates back to Liberman 
in 1977 [7]. The potential due to the OGE appears to be dominated by the contact spin-spin hyperfine interaction (for 
recent reviews see Shimizu et al. [8] and Valcarce et al.[9]). 

Barnes et al. used the "quark Born diagram" formalism to model the dominant contact spin-spin interaction with 
constitute interchange (OGE+CI) to preserve color neutrality. A typical diagram for the process is present in Fig. 2. 
The OGE+CI interaction is an alternative model of short range NN dynamics - compared with the usual heavy meson 
exchange model. Also, it illustrates that the different scales of physics can be self-consistently incorporated into a 
potential using the quark model. 



THE OXFORD POTENTIAL 

For the Oxford potential, we modeled the NN interaction as a combination of one pion exchange (OPE), one sigma 
exchange (OSE) and the OGE+CI potential. For the OPE and OSE, we use the standard the one meson exchange 
mechanism, with the form-factors calculated from the 3 Po model. Charge dependence was incorporated as in CD- 
Bonn^]. 

In order to fit the P and F wave phase shifts, additional repulsion was needed in these channels. Therefore, we added 
one omega exchange in these channels with a pure vector, isovector coupling in these channels. A detailed presentation 
of the potential will be the subject of a forthcoming paper[10] and the quark level origins of the CO exchange effects 
will be investigated. 
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FIGURE 3. Phase shift and inelasticity predictions of the Oxford model for the S, P, and Di partial waves and the E] mixing 
angle. Data from the SAID program[l 1] are the black dots while the predictions are the solid red curves. 

RESULTS 

The parameters of the Oxford potential were fitted to the single-energy values of the NN scattering phaseshifts as 
determined from the SAID analysis program[l 1] and to the properties of the deuteron. No attempt has been made to 
perform a fit to NN scattering data of this 'proof of principle' potential. As the potential does not include Coulomb 
and higher-order electromagnetic corrections, only data in the np scattering channel have been considered in the fit. 



Scattering Phase Shifts 

The NN potential is non-perturbatively strong. Therefore, we solve a partial wave decomposed Lippmann- 
Schwinger type equation. The phase shifts are extracted from the resulting T matrix using standard techniques. We 
present the results of our phase shift parameterization in Fig. 3. The results of the highly accurate potentials are not 
presented - all the potentials reproduce the experimental phase shifts well. 

As can be seen from Fig. 3, the Oxford potential is able to reproduce the low-lying phase shifts very well. We only 
considered data up to 7]^ = 400MeV. Higher L phase shifts are not presented here, but achieve a comparable fit. 



Deuteron Properties 

The deuteron is the only NN bound state. We solve the Schrodinger Equation to calculate the bound state properties 
of the deuteron. The deuteron properties of the Oxford potential are presented along with the results of the other 
highly accurate potentials and the experimental values. Detailed discussion of the experimental values may be found 
elsewhere [4]. The mass of the sigma was fine-tuned to reproduce the binding energy, the other properties are as 
followed without subsequent adjustment. 



Most quantities are accurate to within 1 % when the experimental error is considered. Given the simplicity of the 
model and its constraints, we find these results highly encouraging. In particular we draw attention to the quadrupole 
moment which the Oxford potential is able better reproduce than the highly accurate potentials. 

TABLE 1. The deuteron properties of the Oxford, CD-Bonn, Nijmegen I, and Argonne vl8 potentials 
along with experimental values. Relativistic and Meson Current corrections have not been included. 



Property 


Oxford 


CD-Bonn 


Nijmegen I 


Argonne vl8 


N 3 LO 


Experiment 


B d [MeV] 


2.2246 


2.224575 


2.224575 


2.224575 


2.224575 


2.224575±0.000009 


Aslfm 1 / 2 ] 


0.8918 


0.8846 


0.8841 


.8850 


0.8843 


0.8848±0.0009 


n 


0.0262 


0.0256 


0.0253 


.0256 


0.0256 


0.0256±0.004 


V^ytfm] 


1.9767 


1.966 


1.9666 


1.967 


1.978 


1.971 ±0.006 


QdV™ 2 ] 


0.2871 


0.270 


0.2719 


.270 


0.275 


0.28590±0.00030 


Pd [%] 


5.604 


4.85 


5.664 


5.76 


4.51 





S-wave Scattering Lengths and Effective Ranges 

Our results for scattering lengths, a N , and effective ranges, r™ , due to the strong nuclear force for the S-wave 
channels are summarized in Table 2 in comparison to experiment and predictions of the other highly accurate 
potentials. The channels required the introduction of small charge symmetry breaking in addition to charge 
independence breaking. In the 3 Si channel, the values of the af and rf quoted were obtained using parameters fitted 
only to the binding energy of the deuteron. We find the discrepancy of the order of 1 % between predictions of the 
Oxford potential and experiment encouraging. 

TABLE 2. The Scattering Lengths and Effective Ranges for the Oxford, CD-Bonn, Nijmegen 
I[12], Argonne vl8, and N 3 LO potentials along with experimental values. 

Oxford CD-Bonn Nijmegen I Argonne vl 8 N 3 LO Experiment 



'So 

a N -17.365 -17.4602 -17.164 -17.083 

rf 2.886 2.845 2.865 2.876 

af n -18.949 -18.9680 -18.818 -18.900 -18.9 ±0.4 

rf„ 2.859 2.819 2.834 2.838 2.75 ±0.11 



a N 

y 

of 5.46 5.4196 5.4194 5.402 5.417 5.419 ±0.007 

rf 1.77 1.751 1.7536 1.752 1.752 1.753 ±0.008 



-23.85 -23.7380 -23.084 -23.732 -23.74 ±0.020 

2.753 2.671 2.703 2.752 2.77 ±0.05 

3 S, 



HIGH MOMENTUM BEHAVIOR 

It has been known for some time that the highly accurate potentials remain strong even for high values of momentum 
inducing strong correlations in many-body wave functions. Sophisticated techniques have been developed to integrate 
out the unphysical high momentum degrees of freedom while retaining the physical, low momentum behavior of the 
potentials [13, 14]. 

In Fig. 4 we plot the potential matrix elements in the S and triplet P channels as functions of momentum for 
the Oxford and Nijm I potentials. Because of the quark model wave functions the Oxford potential has a non-local 
Gaussian form-factor which strongly suppresses the potential as momentum increases. This form-factor is similar 
to the one employed in the chiral N 3 LO potential. The rapid fall off in our potential, as opposed to say Nijm I, is 
physically preferred. 



CONCLUSIONS 



We have been conducting research into the possibility of calculating NN interactions mechanisms from the quark 
model. To demonstrate that this technique can lead to an accurate, physically motivated potential, we have produced 
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FIGURE 4. The Oxford and Nijm I on-shell potential matrix elements for a range of momentum. 



the Oxford potential. 

In this paper we reported on the preliminary predictions of the model. In addition to an excellent reproduction 
of the phase shifts, the potential was able to accurately reproduce deuteron and low-energy scattering properties to 
within 1%. Given the simplicity of this model, we find these results to be very encouraging, particularly as regards 
the quadrupole moment of the deuteron. The potential is also the first to exhibit the correct high momentum behavior 
resulting from theoretically motivated form-factors. 

A paper which presents the Oxford potential in detail is under preparation by the authors [10]. Nuclear matter 
calculations have been carried out and show that the Oxford potential has a competitive equation of state. Future 
research into the high momentum behavior of the potential is planned. 
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